Repeated sequences as genetic markers in pooled tissue samples

We show, using the PDR1 element of pea, that dispersed repeated sequences of moderate copy number can be used simply and efficiently to generate markers linked to a trait of interest. Inspection of hybridization patterns of repeated sequences to DNA mixtures of pooled genotypes is a sensitive way of detecting such markers. The large number of bands in tracks of digests of these mixtures allows the simultaneous sampling of loci at many places in the genome, and the many unlinked loci serve as internal controls. It is also shown that intensity ratios calculated from these band differences can be used to give a rough estimate of linkage distance.

Development of type I genetic markers from expressed sequence tags in highly polymorphic species

Expressed sequence tag (EST) databases provide a primary source of nuclear DNA sequences for genetic marker development in non-model organisms. To date, the process has been relatively inefficient for several reasons: - 1) priming site polymorphism in the template leads to inferior or erratic amplification; - 2) introns in the target amplicon are too large and/or numerous to allow effective amplification under standard screening conditions, and; - 3) at least occasionally, a PCR primer straddles an exon–intron junction and is unable to bind to genomic DNA template. The first is only a minor issue for species or strains with low heterozygosity but becomes a significant problem for species with high genomic variation, such as marine organisms with extremely large effective population sizes. Problems arising from unanticipated introns are unavoidable but are most pronounced in intron-rich species, such as vertebrates and lophotrochozoans. We present an approach to marker development in the Pacific oyster Crassostrea gigas, a highly polymorphic and intron-rich species, which minimizes these problems, and should be applicable to other non-model species for which EST databases are available. Placement of PCR primers in the 3′ end of coding sequence and 3′ UTR improved PCR success rate from 51% to 97%. Almost all (37 of 39) markers developed for the Pacific oyster were polymorphic in a small test panel of wild and domesticated oysters.

Novel molecular markers of Chlamydia pecorum genetic diversity in the koala (Phascolarctos cinereus)

Background Chlamydia pecorum is an obligate intracellular bacterium and the causative agent of reproductive and ocular disease in several animal hosts including koalas, sheep, cattle and goats. C. pecorum strains detected in koalas are genetically diverse, raising interesting questions about the origin and transmission of this species within koala hosts. While the ompA gene remains the most widely-used target in C. pecorum typing studies, it is generally recognised that surface protein encoding genes are not suited for phylogenetic analysis and it is becoming increasingly apparent that the ompA gene locus is not congruent with the phylogeny of the C. pecorum genome. Using the recently sequenced C. pecorum genome sequence (E58), we analysed 10 genes, including ompA, to evaluate the use of ompA as a molecular marker in the study of koala C. pecorum genetic diversity. Results Three genes (incA, ORF663, tarP) were found to contain sufficient nucleotide diversity and discriminatory power for detailed analysis and were used, with ompA, to genotype 24 C. pecorum PCR-positive koala samples from four populations. The most robust representation of the phylogeny of these samples was achieved through concatenation of all four gene sequences, enabling the recreation of a "true" phylogenetic signal. OmpA and incA were of limited value as fine-detailed genetic markers as they were unable to confer accurate phylogenetic distinctions between samples. On the other hand, the tarP and ORF663 genes were identified as useful "neutral" and "contingency" markers respectively, to represent the broad evolutionary history and intra-species genetic diversity of koala C. pecorum. Furthermore, the concatenation of ompA, incA and ORF663 sequences highlighted the monophyletic nature of koala C. pecorum infections by demonstrating a single evolutionary trajectory for koala hosts that is distinct from that seen in non-koala hosts. Conclusions While the continued use of ompA as a fine-detailed molecular marker for epidemiological analysis appears justified, the tarP and ORF663 genes also appear to be valuable markers of phylogenetic or biogeographic divisions at the C. pecorum intra-species level. This research has significant implications for future typing studies to understand the phylogeny, genetic diversity, and epidemiology of C. pecorum infections in the koala and other animal species.

Markers of genetic susceptibility in human environmental hygiene and toxicology : the role of selected CYP, NAT and GST genes

Inherited genetic traits co-determine the susceptibility of an individual to a toxic chemical. Special emphasis has been put on individual responses to environmental and industrial carcinogens, but other chronic diseases are of increasing interest. Polymorphisms of relevant xenobiotic metabolising enzymes may be used as toxicological susceptibility markers. A growing number of genes encoding enzymes involved in biotransformation of toxicants and in cellular defence against toxicant-induced damage to the cells has been identified and cloned, leading to increased knowledge of allelic variants of genes and genetic defects that may result in a differential susceptibility toward environmental toxicants. "Low penetrating" polymorphisms in metabolism genes tend to be much more common in the population than allelic variants of "high penetrating" cancer genes, and are therefore of considerable importance from a public health point of view. Positive associations between cancer and CYP1A1 alleles, in particular the *2C I462V allele, were found for tissues following the aerodigestive tract. Again, in most cases, the effect of the variant CYP1A1 allele becomes apparent or clearer in connection with the GSTM1 null allele. The CYP1B1 codon 432 polymorphism (CYP1B1*3) has been identified as a susceptibility factor in smoking-related head-and-neck squameous cell cancer. The impact of this polymorphic variant of CYP1B1 on cancer risk was also reflected by an association with the frequency of somatic mutations of the p53 gene. Combined genotype analysis of CYP1B1 and the glutathione transferases GSTM1 or GSTT1 has also pointed to interactive effects. Of particular interest for the industrial and environmental field is the isozyme CYP2E1. Several genotypes of this isozyme have been characterised which seem to be associated with different levels of expression of enzyme activity. The acetylator status for NAT2 can be determined by genotyping or by phenotyping. In the pathogenesis of human bladder cancer due to occupational exposure to "classical" aromatic amines (benzidine, 4-aminodiphenyl, 1-naphthylamine) acetylation by NAT2 is regarded as a detoxication step. Interestingly, the underlying European findings of a higher susceptibility of slow acetylators towards aromatic amines are in contrast to findings in Chinese workers occupationally exposed to aromatic amines which points to different mechanisms of susceptibility between European and Chinese populations. Regarding human bladder cancer, the hypothesis has been put forward that genetic polymorphism of GSTM1 might be linked with the occurrence of this tumour type. This supports the hypothesis that exposure to PAH might causally be involved in urothelial cancers. The human polymorphic GST catalysing conjugation of halomethanes, dihalomethanes, ethylene oxide and a number of other industrial compounds could be characterised as a class theta enzyme (GSTT1) by means of molecular biology. "Conjugator" and "non-conjugator" phenotypes are coincident with the presence and absence of the GSTT1 gene. There are wide variations in the frequencies of GSTT1 deletion (GSTT1 *0/0) among different ethnicities. Human phenotyping is facilitated by the GST activity towards methyl bromide or ethylene oxide in erythrocytes which is representative of the metabolic GSTT1 competence of the entire organism. Inter-individual variations in xenobiotic metabolism capacities may be due to polymorphisms of the genes coding for the enzymes themselves or of the genes coding for the receptors or transcription factors which regulate the expression of the enzymes. Also, polymorphisms in several regions of genes may cause altered ligand affinity, transactivation activity or expression levels of the receptor subsequently influencing the expression of the downstream target genes. Studies of individual susceptibility to toxicants and gene-environment interaction are now emerging as an important component of molecular epidemiology.

Impacts of traditional husbandry practices on exploitable levels of genetic diversity in cultured 'Tra' catfish (Pangasianodon hypophthalmus) in the Mekong Delta, Vietnam

Sutchi catfish (Pangasianodon hypophthalmus) – known more universally by the Vietnamese name ‘Tra’ is an economically important freshwater fish in the Mekong Delta in Vietnam that constitutes an important food resource. Artificial propagation technology for Tra catfish has only recently been developed along the main branches of the Mekong River where more than 60% of the local human population participate in fishing or aquaculture. Extensive support for catfish culture in general, and that of Tra (P. hypophthalmus) in particular, has been provided by the Vietnamese government to increase both the scale of production and to develop international export markets. In 2006, total Vietnamese catfish exports reached approximately 286,602 metric tons (MT) and were valued at 736.87 $M with a number of large new export destinations being developed. Total value of production from catfish culture has been predicted to increase to approximately USD 1 billion by 2020. While freshwater catfish culture in Vietnam has a promising future, concerns have been raised about long-term quality of fry and the effectiveness of current brood stock management practices, issues that have been largely neglected to date. In this study, four DNA markers (microsatellite loci: CB4, CB7, CB12 and CB13) that were developed specifically for Tra (P. hypophthalmus) in an earlier study were applied to examine the genetic quality of artificially propagated Tra fry in the Mekong Delta in Vietnam. The goals of the study were to assess: (i) how well available levels of genetic variation in Tra brood stock used for artificial propagation in the Mekong Delta of Vietnam (breeders from three private hatcheries and Research Institute of Aquaculture No2 (RIA2) founders) has been conserved; and (ii) whether or not genetic diversity had declined significantly over time in a stock improvement program for Tra catfish at RIA2. A secondary issue addressed was how genetic markers could best be used to assist industry development. DNA was extracted from fins of catfish collected from the two main branches of the Mekong River inf Vietnam, three private hatcheries and samples from the Tra improvement program at RIA2. Study outcomes: i) Genetic diversity estimates for Tra brood stock samples were similar to, and slightly higher than, wild reference samples. In addition, the relative contribution by breeders to fry in commercial private hatcheries strongly suggest that the true Ne is likely to be significantly less than the breeder numbers used; ii) in a stock improvement program for Tra catfish at RIA2, no significant differences were detected in gene frequencies among generations (FST=0.021, P=0.036>0.002 after Bonferroni correction); and only small differences were observed in alleles frequencies among sample populations. To date, genetic markers have not been applied in the Tra catfish industry, but in the current project they were used to evaluate the levels of genetic variation in the Tra catfish selective breeding program at RIA2 and to undertake genetic correlations between genetic marker and trait variation. While no associations were detected using only four loci, they analysis provided training in the practical applications of the use of molecular markers in aquaculture in general, and in Tra culture, in particular.

Association between prostinogen (KLK15) genetic variants and prostate cancer risk and aggressiveness in Australia and a meta-analysis of GWAS data

Background: Kallikrein 15 (KLK15)/Prostinogen is a plausible candidate for prostate cancer susceptibility. Elevated KLK15 expression has been reported in prostate cancer and it has been described as an unfavorable prognostic marker for the disease. Objectives: We performed a comprehensive analysis of association of variants in the KLK15 gene with prostate cancer risk and aggressiveness by genotyping tagSNPs, as well as putative functional SNPs identified by extensive bioinformatics analysis. Methods and Data Sources: Twelve out of 22 SNPs, selected on the basis of linkage disequilibrium pattern, were analyzed in an Australian sample of 1,011 histologically verified prostate cancer cases and 1,405 ethnically matched controls. Replication was sought from two existing genome wide association studies (GWAS): the Cancer Genetic Markers of Susceptibility (CGEMS) project and a UK GWAS study. Results: Two KLK15 SNPs, rs2659053 and rs3745522, showed evidence of association (p, 0.05) but were not present on the GWAS platforms. KLK15 SNP rs2659056 was found to be associated with prostate cancer aggressiveness and showed evidence of association in a replication cohort of 5,051 patients from the UK, Australia, and the CGEMS dataset of US samples. A highly significant association with Gleason score was observed when the data was combined from these three studies with an Odds Ratio (OR) of 0.85 (95% CI = 0.77-0.93; p = 2.7610 24). The rs2659056 SNP is predicted to alter binding of the RORalpha transcription factor, which has a role in the control of cell growth and differentiation and has been suggested to control the metastatic behavior of prostate cancer cells. Conclusions: Our findings suggest a role for KLK15 genetic variation in the etiology of prostate cancer among men of European ancestry, although further studies in very large sample sets are necessary to confirm effect sizes.

Development and characterisation of tri- and tetra-nucleotide polymorphic microsatellite markers for skipjack tuna (Katsuwonus pelamis)

Skipjack tuna (katsuwonus pelamis) (SJT) is the largest tuna fishery in all the major oceans around the world, and the largest marine fishery in Sri Lanka. Knowledge of genetic population structure and effective population size of SJT in the Indian Ocean and other major oceans, however, is still lacking for better management practices and conservation strategies. We developed microsatellite genetic markers using SJT around Sri Lanka in the Indian Ocean, and characterise one tri- and seven tetra-nucleotide microsatellite loci isolated from enriched genomic libraries from SJT, to provide tools for addressing both conservation and fisheries management questions. An analysis of these eight microsatellite markers in two populations of SJT from eastern Sri Lanka (n = 44) and the Maldives Islands (n = 53) showed that all eight microsatellites were polymorphic with an average number of alleles per locus of 11.80 (range 5-27). Expected heterozygosities at marker loci ranged from 0.450 to 0.961. These markers are being used currently to characterise population structure and extent of natural gene flow in SJT populations from the eastern and western Indian Ocean. No significant linkage disequilibrium was detected among any loci pairs.

Genome-wide search for markers associated with bovine spongiform encephalopathy

A genome-wide search for markers associated with BSE incidence was performed by using Transmission-Disequilibrium Tests (TDTs). Significant segregation distortion, i.e., unequal transmission probabilities of alleles within a locus, was found for three marker loci on Chromosomes (Chrs) 5, 10, and 20. Although TDTs are robust to false associations owing to hidden population substructures, it cannot distinguish segregation distortion caused by a true association between a marker and bovine spongiform encephalopathy (BSE) from a population-wide distortion. An interaction test and a segregation distortion analysis in half-sib controls were used to disentangle these two alternative hypotheses. None of the markers showed any significant interaction between allele transmission rates and disease status, and only the marker on Chr 10 showed a significant segregation distortion in control individuals. Nevertheless, the control group may have been a mixture of resistant and susceptible but unchallenged individuals. When new genotypes were generated in the vicinity of these three markers, evidence for an association with BSE was confirmed for the locus on Chr 5.

Genetic studies of common rheumatological diseases

Editorial

Evidence for an X-linked genetic component in familial typical migraine

Migraine is a common complex disorder that shows strong familial aggregation. There is a general increased prevalence of migraine in females compared with males, with recent studies indicating that migraine affects 18% of females compared with 6% of males. This preponderance of females among migraine sufferers coupled with evidence of an increased risk of migraine in first degree relatives of male probands but not in relatives of female probands suggests the possibility of an X-linked dominant gene. We report here the localization of a typical migraine susceptibility locus to the X chromosome. Of three large multigenerational migraine pedigrees two families showed significant excess allele sharing to Xq markers (P = 0.031 and P = 0.012). Overall analysis of data from all three pedigrees gave significant evidence in support of linkage and heterogeneity (HLOD = 3.1). These findings provide conclusive evidence that familial typical migraine is a heterogeneous disorder. We suggest that the localization of a migraine susceptibility locus to the X chromosome could in part explain the increased risk of migraine in relatives of male probands and may be involved in the increased female prevalence of this disorder.

Serum bilirubin concentration is modified by UGT1A1 Haplotypes and influences risk of type-2 diabetes in the Norfolk Island Genetic Isolate

Background Located in the Pacific Ocean between Australia and New Zealand, the unique population isolate of Norfolk Island has been shown to exhibit increased prevalence of metabolic disorders (type-2 diabetes, cardiovascular disease) compared to mainland Australia. We investigated this well-established genetic isolate, utilising its unique genomic structure to increase the ability to detect related genetic markers. A pedigree-based genome-wide association study of 16 routinely collected blood-based clinical traits in 382 Norfolk Island individuals was performed. Results A striking association peak was located at chromosome 2q37.1 for both total bilirubin and direct bilirubin, with 29 SNPs reaching statistical significance (P < 1.84 × 10−7). Strong linkage disequilibrium was observed across a 200 kb region spanning the UDP-glucuronosyltransferase family, including UGT1A1, an enzyme known to metabolise bilirubin. Given the epidemiological literature suggesting negative association between CVD-risk and serum bilirubin we further explored potential associations using stepwise multivariate regression, revealing significant association between direct bilirubin concentration and type-2 diabetes risk. In the Norfolk Island cohort increased direct bilirubin was associated with a 28 % reduction in type-2 diabetes risk (OR: 0.72, 95 % CI: 0.57-0.91, P = 0.005). When adjusted for genotypic effects the overall model was validated, with the adjusted model predicting a 30 % reduction in type-2 diabetes risk with increasing direct bilirubin concentrations (OR: 0.70, 95 % CI: 0.53-0.89, P = 0.0001). Conclusions In summary, a pedigree-based GWAS of blood-based clinical traits in the Norfolk Island population has identified variants within the UDPGT family directly associated with serum bilirubin levels, which is in turn implicated with reduced risk of developing type-2 diabetes within this population.

Prediction of individual genetic risk to prostate cancer using a polygenic score

BACKGROUND Polygenic risk scores comprising established susceptibility variants have shown to be informative classifiers for several complex diseases including prostate cancer. For prostate cancer it is unknown if inclusion of genetic markers that have so far not been associated with prostate cancer risk at a genome-wide significant level will improve disease prediction. METHODS We built polygenic risk scores in a large training set comprising over 25,000 individuals. Initially 65 established prostate cancer susceptibility variants were selected. After LD pruning additional variants were prioritized based on their association with prostate cancer. Six-fold cross validation was performed to assess genetic risk scores and optimize the number of additional variants to be included. The final model was evaluated in an independent study population including 1,370 cases and 1,239 controls. RESULTS The polygenic risk score with 65 established susceptibility variants provided an area under the curve (AUC) of 0.67. Adding an additional 68 novel variants significantly increased the AUC to 0.68 (P = 0.0012) and the net reclassification index with 0.21 (P = 8.5E-08). All novel variants were located in genomic regions established as associated with prostate cancer risk. CONCLUSIONS Inclusion of additional genetic variants from established prostate cancer susceptibility regions improves disease prediction. Prostate 75:1467–1474, 2015. © 2015 Wiley Periodicals, Inc.

Stem cell–related gene expression in clonal populations of mesenchymal stromal cells from bone marrow

Decline in the frequency of potent mesenchymal stem cells (MSCs) has been implicated in ageing and degenerative diseases. Increasing the circulating stem cell population can lead to renewed recruitment of these potent cells at sites of damage. Therefore, identifying the ideal cells for ex vivo expansion will form a major pursuit of clinical applications. This study is a follow-up of previous work that demonstrated the occurrence of fast-growing multipotential cells from the bone marrow samples. To investigate the molecular processes involved in the existence of such varying populations, gene expression studies were performed between fast- and slow-growing clonal populations to identify potential genetic markers associated with stemness using the quantitative real-time polymerase chain reaction comprising a series of 84 genes related to stem cell pathways. A group of 10 genes were commonly overrepresented in the fast-growing stem cell clones. These included genes that encode proteins involved in the maintenance of embryonic and neural stem cell renewal (sex-determining region Y-box 2, notch homolog 1, and delta-like 3), proteins associated with chondrogenesis (aggrecan and collagen 2 A1), growth factors (bone morphogenetic protein 2 and insulin-like growth factor 1), an endodermal organogenesis protein (forkhead box a2), and proteins associated with cell-fate specification (fibroblast growth factor 2 and cell division cycle 2). Expression of diverse differentiation genes in MSC clones suggests that these commonly expressed genes may confer the maintenance of multipotentiality and self-renewal of MSCs.

Development of rapid and highly resolving combinatorial genotyping schemes for Campylobacter jejuni and Campylobacter coli

Campylobacter jejuni followed by Campylobacter coli contribute substantially to the economic and public health burden attributed to food-borne infections in Australia. Genotypic characterisation of isolates has provided new insights into the epidemiology and pathogenesis of C. jejuni and C. coli. However, currently available methods are not conducive to large scale epidemiological investigations that are necessary to elucidate the global epidemiology of these common food-borne pathogens. This research aims to develop high resolution C. jejuni and C. coli genotyping schemes that are convenient for high throughput applications. Real-time PCR and High Resolution Melt (HRM) analysis are fundamental to the genotyping schemes developed in this study and enable rapid, cost effective, interrogation of a range of different polymorphic sites within the Campylobacter genome. While the sources and routes of transmission of campylobacters are unclear, handling and consumption of poultry meat is frequently associated with human campylobacteriosis in Australia. Therefore, chicken derived C. jejuni and C. coli isolates were used to develop and verify the methods described in this study. The first aim of this study describes the application of MLST-SNP (Multi Locus Sequence Typing Single Nucleotide Polymorphisms) + binary typing to 87 chicken C. jejuni isolates using real-time PCR analysis. These typing schemes were developed previously by our research group using isolates from campylobacteriosis patients. This present study showed that SNP + binary typing alone or in combination are effective at detecting epidemiological linkage between chicken derived Campylobacter isolates and enable data comparisons with other MLST based investigations. SNP + binary types obtained from chicken isolates in this study were compared with a previously SNP + binary and MLST typed set of human isolates. Common genotypes between the two collections of isolates were identified and ST-524 represented a clone that could be worth monitoring in the chicken meat industry. In contrast, ST-48, mainly associated with bovine hosts, was abundant in the human isolates. This genotype was, however, absent in the chicken isolates, indicating the role of non-poultry sources in causing human Campylobacter infections. This demonstrates the potential application of SNP + binary typing for epidemiological investigations and source tracing. While MLST SNPs and binary genes comprise the more stable backbone of the Campylobacter genome and are indicative of long term epidemiological linkage of the isolates, the development of a High Resolution Melt (HRM) based curve analysis method to interrogate the hypervariable Campylobacter flagellin encoding gene (flaA) is described in Aim 2 of this study. The flaA gene product appears to be an important pathogenicity determinant of campylobacters and is therefore a popular target for genotyping, especially for short term epidemiological studies such as outbreak investigations. HRM curve analysis based flaA interrogation is a single-step closed-tube method that provides portable data that can be easily shared and accessed. Critical to the development of flaA HRM was the use of flaA specific primers that did not amplify the flaB gene. HRM curve analysis flaA interrogation was successful at discriminating the 47 sequence variants identified within the 87 C. jejuni and 15 C. coli isolates and correlated to the epidemiological background of the isolates. In the combinatorial format, the resolving power of flaA was additive to that of SNP + binary typing and CRISPR (Clustered regularly spaced short Palindromic repeats) HRM and fits the PHRANA (Progressive hierarchical resolving assays using nucleic acids) approach for genotyping. The use of statistical methods to analyse the HRM data enhanced sophistication of the method. Therefore, flaA HRM is a rapid and cost effective alternative to gel- or sequence-based flaA typing schemes. Aim 3 of this study describes the development of a novel bioinformatics driven method to interrogate Campylobacter MLST gene fragments using HRM, and is called ‘SNP Nucleated Minim MLST’ or ‘Minim typing’. The method involves HRM interrogation of MLST fragments that encompass highly informative “Nucleating SNPS” to ensure high resolution. Selection of fragments potentially suited to HRM analysis was conducted in silico using i) “Minimum SNPs” and ii) the new ’HRMtype’ software packages. Species specific sets of six “Nucleating SNPs” and six HRM fragments were identified for both C. jejuni and C. coli to ensure high typeability and resolution relevant to the MLST database. ‘Minim typing’ was tested empirically by typing 15 C. jejuni and five C. coli isolates. The association of clonal complexes (CC) to each isolate by ‘Minim typing’ and SNP + binary typing were used to compare the two MLST interrogation schemes. The CCs linked with each C. jejuni isolate were consistent for both methods. Thus, ‘Minim typing’ is an efficient and cost effective method to interrogate MLST genes. However, it is not expected to be independent, or meet the resolution of, sequence based MLST gene interrogation. ‘Minim typing’ in combination with flaA HRM is envisaged to comprise a highly resolving combinatorial typing scheme developed around the HRM platform and is amenable to automation and multiplexing. The genotyping techniques described in this thesis involve the combinatorial interrogation of differentially evolving genetic markers on the unified real-time PCR and HRM platform. They provide high resolution and are simple, cost effective and ideally suited to rapid and high throughput genotyping for these common food-borne pathogens.